Phase-Sensitive Tetracrystal Pairing-Symmetry Measurements and Broken 
Time-Reversal Symmetry States of High Tc Superconductors 
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A detailed analysis of the symmetric tetracrystal geometry used in phase-sensitive pairing symmetry 
experiments on high Tc superconductors is carried out for both bulk and surface time-reversal 
symmetry-breaking states, such as the d + id' and d + is states. The results depend critically on 
the substrate geometry. In the general case, for the bulk d + id' (or d + is) state, the measured 
flux quantization should in general not be too different from that obtained in the pure d-wave case, 
provided \d'\ -C |d| (or |s| |d|). However, in one particular high symmetry geometry, the d + id' 
state gives results that allow it to be distinguished from the pure d and the d + is states. Results 
are also given for the cases where surface d -I- is or d -I- id' states occur at a [110] surface of a bulk 
d-wave superconductor. Remarkably, in the highest symmetry geometry, a number of the broken 
time-reversal symmetry states discussed above give flux quantization conditions usually associated 
with states not having broken time-reversal symmetry. 

PACS numbers: 74.25.Fy 



This article describes the use of the symmet- 
ric tetracrystal geometry in phase-sensitive pairing- 
symmetry experiments designed to investigate broken 
time-reversal symmetry superconducting states. The ad- 
vantages of the symmetric tetracrystal geometry for such 
experiments will be outlined, and results will be given for 
a number of specific cases. 

All of the hole-doped high temperature superconduc- 
tors (such as YBa2Cu306+2;) measured up to the present 
in phase sensitive pairing symmetry experiments have ex- 
hibited the behavior expected of a superconducting state 
with x'^ — y^ symmetry (see Ref. However, interesting 
questions related to the symmetry of the superconduct- 
ing state of the high Tc superconductors still remain, and 
a number of investigations of possible deviations from 
a pure x^ — y^ symmetry have been and are bjeing un- 
dertaken. For example, it has been suggestecH that at 
low temperatures, an interaction with magnetic impuri- 
ties can cause the d-wave state to become unstable with 
respect to the formation of a d + id' state. (Here d is 
used to denote a state of x'^ — y"^ symmetry, and d' to 
denote a state of xy symmetry.) Also, the anyan mech- 
anism for high Tc superconductivity is knownd to give 
superconducting states of d -f id' symmetry. Attempts 
to interpretau anomalies observedQfl in the thermal con- 
ductivity of certain bismuth superconductors in terms of 
d + id' superconductivity has further heightened interest 
in broken time-reversal symmetry states. Bulk d + is su- 
perconductivity has also been noted as a possibility^. [It 
is known that because the d + id' state (and the d + is 
state) involves order parameters {d and d') with two dif- 
ferent symmetries, the phase transition to a supercon- 
ducting d + id' state can not occur at a single continuous 
phase transition.] Also, the possible existence of surface 
d + is or d + id' states at a [110] surface of a bulk d- 
wave superconductor has been explored theoretically in 



Refs. |[^-13, particularly with reference to the possible 
splitting of the zero bias conductance peak observed in 
NIS (normal to insulator to superconductor) tunneling 
experiments. Experimental evidence for such a splitting 
has been found, at low temperatures, in Rcf. |lj. 

Very recently, evidenceEj of the half integral flux quan- 
tum effect in tricrystal experiments on the electron-doped 
superconductors Nd2-a:Cea;Cu04 and Pr2_a;Cea:Cu04 
has been presented indicating that the superconducting 
state of these materials also has d-wave symmetry. This 
result should be contrasted with one other recent expeE^ 
imental result, namely a thermal conductivity study££l 
on Pr2_2;Cea;Cu04 indicating that there are no propa- 
gating low-energy quasiparticle excitations in this mate- 
rial. These results (in which a material displays both a 
gapless d-wave state and the absence of propagating low 
energy excitations) are extremely unusual, and while per- 
haps not necessarily contradictory, require further study 
if they are to be understood. One possibility is that the 
absence of propagating low energy states implied by the 
thermal conductivity measurements could be accounted 
for if the low energy quasiparticle excitations were local- 
ized by disorder in Pr2-2;Cej;Cu04, as pointed out in Ref. 
|l| (see also Rcf. |l|). 

A question that might be raised about these results 
is to what extent the observations of a half integral flux 
quantum in a tricrystal experiment could, rather than 
implying a pure d-wave state, be consistent with a d + id' 
state (which has a gap and therefore no low energy quasi- 
particle excitations). The reason that this question might 
be asked is that there is at present no detailed analysis 
indicating what the signature of a d -I- id' superconduc- 
tor would be in a tricrystal experiment. Furthermore, as 
noted above, theoretical arguments have been made that 
the broken time-reversal symmetry d -I- id' state could 
play a role in high Tc superconductivity. Clearly some 
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analysis of the signature of d + id! superconductivity in 
phase-sensitive pairing symmetry experiments would be 
useful in resolving these questions. 

The above discussion indicates considerable current in- 
terest in, and motivation for, the study of various types 
of broken time-reversal symmetry states. The main ob- 
jective of this article is therefore to study how broken 
time-reversal symmetry superconducting states of d + id' 
and d+is symmetry, in both their bulk and surface forms, 
would manifest themselves in phase-sensitive fiairing- 
symmetry experiments of the tetracrystal typeu. (For 
reasons discussed in the following paragraph, the focus 
will be on the tetracrystal geometry only, and no discus- 
sion will be given of the tricrystal geometry, which is the 
geometry that has been used in most of the phase sensi- 
tive pairing symmetry experiments carried out to dateEJ.) 
Some aspects ijeiated to the bulk d+is state have already 
been discussedta. The details of the d+id' case are some- 
what different, however. As a result, the d + id' and the 
d + is states give quite different results for an appropri- 
ately chosen experimental geometry and these differences 
provide a method of experimentally identifying a state of 
d + id' symmetry. 

The symmetric tetracrystal configuration proposed in 
Ref. |l^ and realized experimentally in Refs. 1^,^ is 
shown in Fig. 0. This is called a symmetric geometry be- 
cause the X axis of the figure is an axis of reflection sym- 
metry, and this symmetry plays an essential role both 
in simplifying the analysis, and in allowing an unam- 
biguous interpretation of the experimental results. The 
reason for choosing to analyze the tetracrystal geome- 
try (rather than the tricrystal geometry) in this article is 
that this geometry allows a determination to be made of 
order parameter symmetry from the observed flux quan- 
tization (say half integral) based on symmetry arguments 
only. In the tricrystal experiment, symmetry arguments 
can not be used because the experimental arrangement 
has no useful symmetry elements (such as a plane of re- 
flection symmetry). It is therefore necessary to appeal 
to a phenomenological model for the dependence of the 
Josephson critical current on the misorientation angles 
Qi and Q2 of the two grains on either side of the grain 
boundary Josephson junction. For example, an appropri- 
ate model accounting for the known experimental results 
on cuprate superconductors would be 



(1) 



Here C(Oi,02) is a positive factor taking account of 
the changes in grain boundaey, microstructure with in- 
creasing misorientation anglecil; this effect is expected 
to dominate the overall magnitude of the Josephson 
critical current, which is found experimentally to vary 
roughly as exp[— (9i -|-j^2)/0o] where 89 « 5° for hole- 
doped superconductorsES and Oo ~ 2° for electron-doped 
superconductors^. The factor fsym can change sign 
and reflects the pairing symmetry of the superconduct- 
ing state. For the case of d-wave symmetry the simplest 



choices are fsym = cos (2ei)cos(2e2)0 apd the maxi- 
mum disorder limit fsym — cos(29i H- 262)c3. A system- 
atic study of a variety of different tricrystal geometries 
has established the reliability of this approach (see Ref. 
|l| for a review). In addition, the results of the tricrystal 
experiments are in agreement with those obtained us- 
ing the symmetric tetracrystal geometryld^Ot3 which, as 
noted above, does not refer to a detailed model of the 
angular dependence of the Josephson critical current. 

The exploration of various broken time-reversal sym- 
metry states carried out in this article is more complex 
than testing for d-wave versus s-wave symmetry, since 
both the d + id' and d + is states involve linear com- 
binations of order parameters of two different symme- 
tries. Determining a quantitative formula for the an- 
gular dependence of the Josephson critical current suf- 
ficiently accurate for use with the tricrystal geometry 
is thus even more difficult than for the purely d-wave 
and s-wave cases. The symmetric tetracrystal geometry 
presents no such difficulties, as accurate predictions can 
be made on the basis of symmetry arguments only pro- 
vided the cuprate under study is tetragonal. 
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FIG. 1. The tetracrystal geometry The substrate is 

made from two tetragonal crystals of the same material (e.g. 
SrTiOa), SI and S2, with their c-axes normal to the plane of 
the paper, and their tetragonal a and b axes oriented at 45 
degrees to each other, as indicated by the cross ruling in the 
figure. For technical reasons, the configuration in the figure 
is easiest to manufacture if the crystals SI and S2 are each 
composed of two separate crystals that are fused together 
along the line that is the x-axis in the figure. For this reason, 
this is called the tetracrystal geometry. A superconducting 
layer is deposited epitaxially on the substrate, and the grain 
boundary A'AVBB' forms a Josephson junction. The dashed 
circle is not part of the experimental setup, but its arcs give 
the paths of a contour intergrals described in the text. 
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In fact, for the special case of a = 45° in the symmet- 
ric tetracrystal geometry, an additional symmetry comes 
into play resulting in predictions of remarkable simplic- 
ity (see Table | below) . Thus the symmetric tetracrystal 
geometry, which has certain advantages relative to the 
tricrystal geometry even in the casas-|Wfeere no broken 
time-reversal symmetry is envisagecHlljH, should be the 
geometry of choice for the investigation of differences be- 
tween broken-time reversal symmetry states. 

The nature of the spontaneous flux generated at the 
vertex V in Fig. |^ can be determined by using Ginzburg- 



18 2g,2l The su- 



Landau theory, as described in Refs. 
perconducting order is characterized by an order param- 
eter having two components ipd and ipd' ■ (This is for the 
d + id' case; for the d + is case the two components are 
called i/^d and V's.) Also, the value of the A'th component 
(A is d or d') of the order parameter at the point J (J is 
A or B) in the j'th superconductor (j is 1 or 2 for SI or 

S2) is denoted by V'ij • 

The free energy per unit length of the grain boundary 
Josephson junction at point J is 



AA' J 



A, A' 



V'ij*V'vjea;p ii^ 



J2 



A-d£ + c.c. 



,11 



(2) 



which is consistent with gauge invariance and time- 
reversal symmetry. The line integral from point Jl 
to J2 in this equation is across the grain boundary at 
point J from superconductor SI to S2. The fact that 
the X-axis is an axis of reflection symmetry imposes 
the conditions C\\a = —C\\b and Caam = Cwb for 
A ^ A'. (For the d + is case the corresponding relations 
are CddA — ~CddB,CssA = CssB,CsdA = CgdB, and 
CdsA — —CdsB-) It is these symmetry relations that al- 
lows a definitive interpretation of experiments performed 
in the tetracrystal geometry to be made without recourse 
to a detailed model of the_angular dependence of the 
Josephson junction currents. 

The assumption of a d -f- id' state is implemented by 
taking the order parameters to have the form 



(3) 



Note the additional factor of i (which characterises the 

d + id! state) in the definition of V'^'^ j- Also, ajj can 
have the values ±1 corresponding to the two possible 
variants d±id' , each of which has the same free energy. 
For a bulk broken time-reversal symmetry state in, say, 
superconductor SI, this state would be expected to be 
characterized by either the order parameter d -\- id' or 
d — id' throughout the entire crystal SI because a change 
from one variant to the other would cost energy. Thus 
we would have aAi = o'bi, i-e- that same state of super- 
conductor SI at both junctions. 

On the other hand, it is also interesting to consider the 
case of a superconductor that has pure d-wave symmetry 



in the bulk, but in which the d-wave superconductivity is 
suppressed at the junction interface with the simultane- 
ous appearance of a d' order parameter there, so that at 
the junction interface one has d ± id' superconductivity. 
In this case one might have d+id' surface superconductiv- 
ity in superconductor S2 at junction A and d — id' surface 
superconductivity in superconductor S2 at junction B, so 
that aA2 and aB2 are not necessarily equal. 

Now put the above expressions for the order parame- 
ter into the junction free energy of Eq. |2[ and define a 
critical current per unit length along the grain boundary, 
I A, and a phase 9 a, characteristic of the grain boundary 
Josephson junction at A by 



Attc 



Iab 



9dd + AlC^ A2gd'd' + A2gdd' - <^Aigd', 



where 



.9AA' 



(4) 



(5) 



Similarly, for point B on the grain boundary, define Ib 
and Ob 



47rc 



I Be 



iBE 



gdd + CTBi(TB2ffrf'd' 



i{<^B2gdd' - CTBigd'd)- 

(6) 



The fact that the x axis of Fig. |i| is a line of reflection 
symmetry has the important consequence that the same 
coefficients gw that appear in Eq. ^ also appear in Eq. 
^. The grain boundary junction free energies per unit 
length at points A and B can now be written 



$0 

Fa = - — lACos{(f>A + Oa), 

ZTTC 



and 



where 



Fe 



27rc 



'.72 



Ibcos{4>b +9b +Tr)- 



J2 



— / A-de. 



Jl 



2tt 

$0 



(7) 



(8) 



(9) 



Note that the values of ajj in Eqs. |^ and |^ have not yet 
been chosen. These should be chosen so that the critical 
currents I a and Ib have their maximum possible magni- 
tudes, thus leading ultimately to the lowest possible free 
energy for the system. 

A measurement of the spontaneous magnetic flux 
through the grain-boundary Josephson junction of Fig. 
|l| in a direction normal to the page shows that the flux 
is concentrated in the neighborhood of the vertex V, and 
falls exponentially to zero as one proceeds along the grain 
boundary in the direction of either A or BEI. The charac- 
teristic length scale of this fall off is the Josephson pene- 
tration depth. It is assumed that the points A and B are 
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many Josephson penetration depths away from the ver- 
tex V so that there is no fLux through the grain boundary 
at points A and B, nor is there any current flowing across 
the grain boundary at these points. The free energies Fa 
and Fb will therefore have their minimum values at these 
points, leading to the results 

(l)A + eA + TT ^0, (I)b + Ob=0, {mod2n). (10) 

Next note that since the dashed circle is many penetra- 
tion depths from the vertex V, there are no currents flow- 
ing in the superconductor in its vicinity. It follows that 



27r 



A] 



JB] 



(11) 



where the line integral follows an arc of the dashed circle 
in Fig. |l|. Finally, combining Eqs. ^ |l^, and |ll| leads to 
the result that the spontaneous flux through the grain 
boundary junction at the vertex V is 



spont 



2tt 



(12) 



For the sake of definiteness, choose 6a and 9b (which 
previously have been defined mod 2tt) such that < 
{9b — 9a) < 2tt. Also, up to now the magnetic energy 
i'^/(2L) has been assumed to be small in comparison 
with the free energies of Eqs. ^ and ||, and has been left 
out of the problem. Since the magnetic free energy is a 
positive definite function of the flux <&, the state (i.e. the 
value of n in Eq. |l^) giving the lowest value of the total 
free energy is the one for which the flux is 



2tt 



$0- 



(13) 



This minimum flux is such that — $o/2 < $min < $o/2- 
Thus the state corresponding to n = in Eq. ^ is the 
ground state, and the states corresponding to n a nonzero 
integer are possible metastable states. 

In discussing the bulk d ± id' states, the requirements 
aAj — cTBj (see above) are inserted into Eqs. ^ and ||, 
yielding 9b — —9a- An interesting limit is that in which 
the d component of the order parameter is assumed to 
be much greater than the d' component. An examina- 
tion of Eqs. H and |^ in this case shows that, so long as 
a is not too close to 45°, the term in gdd is expected to 
dominate on the right hand side, and hence \9a\ — 
will be much less than tt. This results in an approx- 
imately half integral flux quantum effect, which is not 
very different from what is expected for a pure d state. 
The departure of the spontaneously generated flux from 
$0/2 will grow as the magnitude of -^d' grows relative to 
that of V'd- Thus, to detect broken time-reversal symme- 
try with iV'd'l I'f/'dl requires a high precision measure- 
ment of the spontaneous flux, because the signature of 
the presence of a small d' component of the order param- 
eter is a small departure of the spontaneous flux from 



$0/2. For a = 45° however, an additional symmetry of 
the Josephson junctions requires that gdd = gd'd' = 0. 
This yields 9a = —9b = ±7r/2 and a spontaneous flux 

= n^o, which is precisely the same result as would be 
obtained for a classical s-wave superconductor. The re- 
sult that the broken time-reversal symmetry d + id' state 
has $ = n^o as its flux quantization condition for this 
geometry is quite remarkable since broken time-reversal 
symmetry states are generally believed to have the flux 
quantization condition $ = (n + f)^o, where the fraction 
/ satisfies —1/2 < / < 1/2 and is neither precisely nor 
precisely 1/2. In summary note that, for the d+id' state, 
the geometries with a ^ 22° and a = 45° lead to com- 
pletely different flux quantization conditions, and that 
these different conditions give a way of experimentally 
identifying a d + id' state. 

The introduction to this article cites several theoret- 
ical studies that suggest a bulk rf-wave superconductor 
will tend to exhibit a suppressed c?-wave order parameter 
at a [110] surface, possibly accompanied by the presence 
of a d' or s component at this surface, and with these 
components phased in such a way that the overall state 
at the surface has d-\- id' or d-\-is symmetry. An appro- 
priate geometry in which to study such states is to take 
the angle a in Fig. |^ to be 45 degrees. There is then the 
possibility of a surface d+id' or surface d+is state in the 
superconductor S2 (but not SI) along the grain bound- 
aries A'AV and VBB'. There are two inequivalcnt possi- 
bilities for the surface states of the two grain boundary 
segments: A'AV can be d±id' while VBB' is also d±id', 
or A'AV can be d ± id' while VBB' is d =F id' (here take 
only the upper signs, or only the lower signs). There 
is an interaction energy at the vertex V which should 
give these two inequivalcnt possibilities different free en- 
ergies, but this interaction energy will be much smaller 
than the Josephson energy of the entire grain boundary 
and therefore the assumption will be made below that 
both possibilities can occur for the same sample. 

Prom the work already done above, it follows that the 
case of purely d-wave symmetry in the bulk, with d ± id' 
surface superconductivity at the [110] surface of super- 
conductor S2 for the case of a = 45° also yields a flux 
quantization condition $ — n$o for the case aA2 — o'b2, 
but with the new feature that there is the possibility of 
a states of the same sample having the flux quantization 
condition $ = (n-|- ^)$o for (7^2 = —o'B2- Of all of these 
states, the ground state is the one with zero flux. Again, 
it is remarkable that the broken time-reversal symmetry 
surface states for this geometry lead to flux quantization 
conditions usually thought to be characteristic of states 
without broken time reversal symmetry. 

The analysis of the d + is case is similar, except that 
the deflning equations for the critical currents and char- 
acteristic phases, Eqs. ^ and ^, are replaced by 



Attc 



IaS^^^ = gdd - CTAlO-A2ffss + ii-(^Aigsd + CrA2gds 



(14) 
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and 



+ cfBiOB2gs8 + i{crBigsd + crB2gds)- 



(15) 

Arguments similar to those given above lead to the result 
that the spontaneous flux at the vertex V in Fig. |l| is 
again given by Eq. ^ 

Now consider the limit of the bulk diis case in which 
the magnitude of the d component of the order parameter 
is much greater than that of the s component. An exam- 
ination of Eqs. |l^ and |l^ in this case shows that, so long 
as a is not too close to 45°, the term in g^d is expected 
to dominate the right hand side, and hence \6a\ and \9b\ 
will be much less than tt. This results in an approx- 
imately half integral flux quantum effect, which is not 
very different from what is expected for a pure d state, 
nor from the d + id' case for the corresponding values of 
a. Now note that for a = 45°, the additional symmetry 
of the Josephson junctions requires that gdd = gsd = 0. 
(Also, \gss\ < \gss\-) This yields Oa ^ Og 7r/2 or 
6a — Ob — —tt/2. Hence for a — 45° also the approxi- 
mate half integral flux quantum effect is preserved, which 
is very different from what was predicted above for the 
bulk d + id' case. Thus, the symmetric tetracrystal ge- 
ometry for a = 45° can in principle distinguish between 
the bulk d + id' and d + is states. 

The case of purely d-wave symmetry in the bulk, with 
d ± is surface superconductivity at the [110] surface of 
superconductor S2 for the case of a = 45° yields the flux 
quantization conditions of $ = n^Q for 17^2 = crB2 and 
$ = (n + ^)<i>o for aA2 = — o'_B2- These and other results 
are summarized in Table ^ 

The grain boundaries that occur in Fig. for a = 45° 
are called 45° [001] tilt grain boundaries. Scanning 
SQUID microscope imaging of such grain boundaries 
in YBa2Cu306-i-a; reveals very weak spontaneously gf 
erated delocalized flux along these grain boundaries 
This is thought to be due to either grain-boundary 

TABLE I. Tfie flux quantization conditions for the 
tetracrystal geometry of Fig. ^ with the angle a = 45° for 
a number of different broken time-reversal symmetry states. 
For the surface case, the two states given represent the sur- 
face states in superconductor S2 on grain boundaries A and 
B, respectively, of Fig. Also, / satisfies — ^ < / < |, but is 
not zero; although / can not be precisely i, it is expected to 
be close to i if the time-reversal symmetry breaking is weak. 



state {a — 45° in all cases) 



flux quantization 



bulk d + id' 
bulk d + is 

surface d + id' , d + id' 
surface d + id' , d~ id' 
surface d + is, d + is 
surface d + is, d — is 



n$o 

(n + /)<E>o 

(n-hi)$o 
(n-hi)$o 



facetingE3, or the existence of orthorhombic twinsB^ in the 
YBa2Cu306-i-a;, or both. To avoid this effect (which can 
be detected using scanning SQUID microscopy as just 
mentioned) the material studied should be tetragonal and 
not orthorhombic (thus avoiding twinning) and the grain 
boundaries should be as straight as possible. In the case 
of d+is and d+id' states, the additional component of the 
order parameter (in addition to the d-wave component) 
provides an additional coupling which, if large enough, 
will overcome the effects giving rise to the spontaneous 
delocalized grain boundary flux (which results from the 
fact that usual c?-wave Josephson coupling is zero for an 
ideally flat 45° [001] tilt grain boundary). 

This article has considered measurements of the pair- 
ing symmetry made using the tetracrystal geometry. For 
the bulk d+id' and d+is states, and for the value a ^ 22° 
(see Fig. |l|) advocated in Refs. ^ and |l^, the flux quan- 
tization condition is $ — (n + f)^o, where / is a fraction 
satisfying —1/2 < / < 1/2; also / ^ 0. If the magnitude 
of the symmetry breaking is small, / is expected to be 
close to 1/2, and the measurements must be relatively 
high precision to detect the small differences of the spon- 
taneous flux from the value of $o/2. However, striking 
differences between the bulk d+id' and d+is states occur 
in the highest symmetry case where the angle a = 45°. 
In this case, the d+id' state should exhibit a spontaneous 
flux of zero while the d + is state should exhibit a spon- 
taneous flux of approximately ^o/2. This should provide 
a way of experimentally identifying a state with d + id' 
symmetry. For bulk superconductivity that is purely d- 
wave, and surface superconductivity of the type d+id' or 
d + is at junctions A and B in superconductor S2 of Fig. 
|l| with a = 45° , the flux quantization may be either of 
the type $ = n<i>o or of the type $ = (n-l- ^)<i>o, depend- 
ing on whether the surface states on the two differently 
oriented grain boundaries in Fig. |^ are the same, or are 
time-reversed conjugates. Interestingly, in this highest 
symmetry geometry, flux quantization conditions usually 
expected for states without broken time-reversal symme- 
try [i.e. = $0 or $ = (n -I- 5)^*0] are also found for 
several broken time-reversal symmetry states. 

I would hke to thank J. Y. T. Wei and C. C. Tsuei for 
the discussions that led to this article, C. C. Tsuei and 
L. Taillefer for useful comments on the manuscript, and 
the Natural Sciences and Engineering Research Council 
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